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ABSTRACT 

The dynamics of stars in the inner regions of nearby globular clusters (GCs) 
such as Gl indicate the presence of central concentrated dark masses, and one 
would like to know whether these are indeed intermediate mass black holes 
(IMBHs). As the number of surrounding stars, and their motions, are roughly 
known, the capture rate can be estimated; the question then arises of whether 
the apparent quiescence of the nuclei of these GCs is compatible with the IMBH's 
presence. The role of debris from disrupted stars in activating quiescent nuclei of 
GCs is studied here employing three-dimensional hydrodynamics simulations. It 
is argued that when individual stars are disrupted, the bulk of the debris would 
be swallowed or expelled rapidly compared with the interval between successive 
disruptions. A transient (predominantly of soft X-ray emission) signal could per- 
sist steadily with L ~ i^Edd = 10 41 (M h /10 3 M o ) erg/s for at most tens of years; 
thereafter the flare would rapidly fade. While the infall rate declines as t~ 5 ^ 3 , 
some material may be stored for a longer time in an accretion disk. The IMBH 
luminosity could then remain as high as Lx > 10 39 erg/s for several hundreds of 
years after disruption. In a given object, this ultraluminous X-ray activity would 
have a duty cycle of order 1CT 4 . Quiescent GCs, those with IMBHs now starved 
of fuel, should greatly outnumber active ones; observational constraints would 
not then be stringent until we had observed enough candidates to constitute a 
proper ensemble average. 

Subject headings: black hole physics - hydrodynamics - globular clusters: general 

1. Introduction 

Suggestive evidence has accumulated that intermediate mass black holes exist in some 
globular clusters. There is some dynamical evidence for a mass concentration within the 
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central regions of some globular clusters. The dynamics of stars in the inner regions of 
nearby clusters such as M15, Gl and u Centauri suggest the presence of black holes with 



masses of about 10 3 , 10 4 and 4 x 10 4 M Q , respectively ( 


Gerssen et al. 


2002, 


2003; 


Gebhardt 


et al. 


2002 


2005; 


Pooley &; Rappaport 


2006 


Noyola et al. 


2008 


). The peculiarities of an 



unresolved radio source in Gl indicate some unique object at the center (Ulvestad et al. 



2007). Recently, somewhat arguable evidence has arisen for the presence of IMBHs in young 



star clusters, where ultraluminous, compact X-ray sources (ULXs) have been preferentially 



found to occur (Zezas et al. 2002). Their high luminosities suggest that they are IMBHs 
rather than binaries containing a normal stellar mass black hole ( Portegies Zwart et al.|2004 ). 
But before accepting this conclusion (and dismissing alternative ideas) one would like some 
independent corroboration of the IMBH hypothesis, or, conversely, some way of ruling it 
out. 

We are used to the idea that black holes are implicated in the most powerful sources in 
the universe, and can (when accreting) be ultra-efficient radiators. But there is, for example, 
no sign of such activity in M15 - the X-ray upper limit is no more than 6 x 10 32 erg/s (Ho 



et al. 2003). So could a black hole be so completely starved of fuel that it does not reveal 



its presence? We do not directly know how much gas there is near the IMBH (Pfahl & 



Rappapo"rt||2001 Pooley & Rappaport 2006), and there is no a priori reason why this region 



should be swept clean of gas. Gas that is lost from nearby stars (Baumgardt et al. 2006) 



as well as mass transfer from potentially bound companions (Blecha et al. 2006; Patruno et 



al. 2006 Hopman et al. 2004) can produce observable signatures. However, it is uncertain 
whether gas can be adequately supplied to explain ULX activity. The star density, on the 
other hand, is much better known - after all, if the stars were not closely packed near the 
center of the cluster we would not have evidence for the black hole at all. As stellar orbits 
diffuse in phase space, it therefore seems inevitable that some may wander sufficiently close 
to the hole that they suffer tidal disruption. When a star is disrupted, there is bound to be 
some radiation from the sudden release of gas. The flares resulting from a disrupted star 
could be the clearest diagnostic of a IMBH's presence. 

It is an intricate although tractable problem in stellar dynamics to calculate the chance 



that a star passes within the IMBH's tidal radius (Frank & Rees 1976). In a simple case 



when the velocities are isotropic, the frequency with which a star would enter the zone of 
vulnerability is £ ~ 10™ 7 M^3 3 (n*/10 6 pc~ 3 )(a/10 km s -1 ) -1 yr _1 , where is the stellar 
number density in the cluster nucleus and M^g is the black hole's mass in units of 1O 3 M . 
This estimate although simplified, agrees well with the fiducial rates derived from detailed 
N— body simulations of multi-mass star clusters containing IMBHs ( Baumgardt et aL||2004 ). 
It may seem, however, that even the modest rate of stellar disruptions given above could 
have conspicuous consequences. The debris from a disrupted one solar-mass star per hundred 
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million years, swallowed steadily with 10 per cent radiative efficiency, would yield a luminos- 
ity of L s teady ~ 6 x 10 37 erg/s - higher than is observed for X-ray binaries in quiescence. But 
in reality, one expects brighter flares with short duty cycles, as the time it takes to digest or 
expel the debris from one star is much shorter than the mean interval between one stellar 
disruption and the next. This Letter is concerned with the observational manifestations of 
such phenomena, with particular reference to globular clusters where the masses of the black 
holes (if indeed present) are perhaps of order 10 3 — 1O 4 M . 



2. The tidal Disruption of solar-type star by an IMBH 



A star interacting with a massive black hole cannot be treated as a point mass if it 
gets so close to the hole that it becomes vulnerable to tidal distortions. Such effects become 
important when the pericenter distance becomes as small as the tidal radius: R? ~ 5 x 
lO 11 M h V 3 3 ( J R Jl ,/ J R )(M,/M )- 1 / 3 cm. The gravitational radius, R g ~ 1.5 x 10 8 M hj3 cm, scales 
with mass, whereas the tidal radius goes only as the cube root. The tidal forces at their 
"surfaces" are thus more gentle for black holes of larger mass, and solar-type stars would 
be disrupted only after passing irreversibly inside an ultra-massive hole's horizon: < 
7 x 10 7 M . 

When a rapidly changing tidal force starts to compete with a star's self-gravity, the 
material of the star responds in a complicated way, being stretched along the orbital direc- 



tion and squeezed at right angles to the orbit ( Carter & Luminet 1983 ; Rees 1988 ; Evans & 



Kochanek 1989). To study this problem, we use a three dimensional smoothed particle hy- 



drodynamics method (SPH) to solve the equations of hydrodynamics. Due to its Lagrangian 
nature SPH is perfectly suited to follow tidal disruption processes during which the corre- 
sponding geometry, densities and time scales are changing violently. The SPH-formulation 
that we use in this study is described in Rosswog et al. ( 2008af b). Common to all runs is 
the initiation of the calculations with the star being place safely outside Rt and set it onto 
a parabolic orbit so that -R m i n = Rt/3 for a 1O 3 M black hole. We have considered three 
different initial conditions for the approaching star, constructed here by solving the spher- 
ically symmetric Lane-Emden equations: A [1, 1, 0.6] M solar-type star modeled with a 
polytropic equation of state with adiabatic index T = [5/3, 1.4, 1.4] and i?* = [1, 1, O.75]i? . 

The gross quantitative behavior of a solar-type star plunging deeply within the tidal 
radius is examined here (the reader is refer to ^3] for a discussion on the role of the stellar 
density structure in shaping the history of the mass accretion rate). Several snapshots 
taken from our numerical simulations of a 1 M solar-type star (modeled with T = 5/3) 
are shown in Figure [I] The tidal bulge raised on the star by the black hole becomes of 
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an order unity distortion near pericenter. The resultant gravitational torque spins it up 
to a good fraction of its corotation angular velocity by the time it gets disrupted. This 
takes place on a timescale comparable to the crossing time of the star through periastron, 
At ~ Rt/v p ~ 349.6(M,/M )- 1 / 6 ( J R,/ J R ) 3 / 2 M~3 1/3 s. 

The energy required to tear the star apart (that is the star's self binding energy) is of 
order M*t> 2 , where t>* = (GM*/ R*) 1 ' 2 . During tidal disruption, this energy is supplied at the 
expense of the orbital kinetic energy, which at pericenter ~ Rt is larger by ~ (Mh/M*) 2 / 3 



(Rees 1988). Figure |2j displays the evolution of the differential mass distribution in specific 
energy for the debris. Disruption reduces the orbital energy by the binding energy of the 
star e* ~ (GM+/R*) ~ 10~ 5 c 2 , which is much smaller than the specific kinetic energy at 
pericenter. The variation of the specific energy in the released gas is determined mainly 
by the relative depth of a mass element across the disrupted star in the potential well of 

1 /3 

the black hole. The spread in this specific energy is of order vAv ~ 10 M h3 c 2 , where 
v m v p = c[R g / Rt} 1 / 2 and Av fa v* = [GMj/fl,] 1 / 2 . This is much larger than e* and, as a 
result, almost half the debris escapes on hyperbolic orbits with speeds ~ 3000 M^fkms -1 ; 
the kinetic energy output being ~ 10 50 Mj| g 3 erg (comparable to the energy of a supernova). 
The material would be concentrated in a fan close to the orbital plane. Adiabatic cooling, as 
the material expands, severely reduces the internal radiative energy content before the debris 
became optically thin. The escaping radiation from the unbound debris would therefore 
release much less than the initial energy content of the star. There would therefore be no 
conspicuous flare, until the bound debris fell back onto the IMBH. 



3. Ingesting the Stellar Debris 

The returning gas does not immediately produce a flare of activity from the black 



hole. First material must enter quasi-circular orbits and form an accretion torus (Evans 



& Kochanek 1989). Only then will viscous effects release enough binding energy to power 
a flare. The bound orbits are very eccentric, and the range of orbital periods is large. The or- 
bital semi-major axis of the most tightly bound debris is a ~ 1O 4 M- 3 1/3 ( J R,/ J R )(M,/M )- 2 / 3 R s , 

and the period is only t a ~ 6300(a/1.5 x WR g fl 2 M~l /2 s. If the gaseous debris suffered 
no internal dissipation due to high viscosity or shocks, it would, after one or two orbital 
periods, form a highly elliptical disk with a big spread in apocentric distances between the 
most and least bound orbits, but where at pericenter the orbits are all squeezed in a range 
SR/R ~ R*/R T = O.^Mfc/Mh^) 1 / 3 . As the stream approaches pericenter, the radial focusing 
of the orbits therefore acts as an effective nozzle (Figure |3|. After pericenter passage, the 
outflowing gas is on orbits which collide with the infalling stream near the original orbital 
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plane at apocenter, giving rise to an angular momentum redistributing shock (Figure [3]) 
much like those in cataclysmic variable systems. The debris raining down would, after little 
more than its free-fall time, settle into a disk. 

This orbiting debris starts forming a disk when the most highly bound debris falls back. 
The simulation shows that the first material returns at a time < t a , with a peak infall rate 
roughly given by 31M yr _1 (Figure [I]). Such high infall rates are expected to persist, 
relative steadily, for at least a few orbital periods, before all the highly bound material rains 
down. Since the amount of highly bound material [dMj de) depends sensitively on the density 
structure of the star, we find that more centrally condensed stars produce mass feeding rates 
that are slightly larger in comparison to less centrally condensed stars of comparable mass 
although they reach the self-similar oc t~ 5 / 3 phase after fewer orbital periods. The vicinity 
of the hole would thereafter be fed solely by injection of the inf ailing matter at a ratd^that 
drops off roughly as t _5//3 for t > t&> ~ 10 5 s. Once the torus is formed, it will evolve under 
the influence of viscosity, radiative cooling winds and time dependent mass inflow. 

The black hole cannot accept matter, with a radiative efficiency e = O.leo.i, at a rate 
exceeding MEdd ~ 2 x 10~ 5 eo.iMh,3M Q yr~\ without exceeding the black hole's Eddington 
luminosity. The rate at which the stellar debris returns to the vicinity of the black hole 
exceeds this limit: M ~ 10 6 MEdd for £ = 0.1. Given that this high rate can only be 
sustained by infall of the highly bound material for a time tfb, we infer that, if the viscosity 
were high enough to process all the material within the infall timescale, then the luminosity 
of the hole could not remain as high as the Eddington luminosity for longer than ^Edd ~ 
5 x 10 3 e ^ 5 M h g^ 5 to, ~ 17e Q ^ 5 M h l yr. The viscosity would have to be implausibly low 
(i.e. the usual viscous dissipation time t& for a thick disk would be a -1 ; a would have to 
be below 1CT 5 for t d > i E dd) f° r the bulk of the mass to be stored for longer than £ E dd 
in a reservoir at R ~ R^- A luminosity ~ i^Edd — 10 41 Mh 5 3 erg s _1 can therefore only be 
maintained for at most ten years; thereafter the flare would continue to fade as t~ 3 ^ 5 . It 
is clear from the behavior of M that most of the debris would be fed to the hole far more 
rapidly than it could be accepted if the radiative efficiency were high; much of the bound 
debris must either escape in a radiatively-driven outflow or be swallowed inefficiently 



1 Note that the rate at which mass returns to the black hole is only proportional to t~ 5 ^ 3 when dM/de is 
constant (Rees 1988 Phinney~| |l989). Figure [2] displays the differential mass distributions in specific energy 
obtained by the simulation, which show that dM/de is nearly constant. 



4. Discussion 



4.1. Observability of L Edd Flares 

A distinctive consequence of a 10 3 — 10 4 M & IMBH's presence in the centers of globu- 
lar clusters would be the transient flares produced as the bound debris from the disrupted 
solar-type stars is swallowed, the luminosities being as high as L Edd = 10 41 M h 3 erg s -1 . The 
rise and the peak bolometric luminosity can be predicted with some confidence. However, 
the effective surface temperature (and thus the fraction of luminosity that emerges predom- 
inantly in the soft X-ray band) is harder to predict, as it depends on the size of the effective 
photosphere that shrouds the hole. For t < £ Edd , the effective surface temperature should be 
T < T Edd = (L Edd /47rcr_Rg) 1 / 4 ~ 1M^ A keV. In a given globular cluster, these flares would 
have a duty cycle of order £t Edd ~ 10~ 6 (£/10~ 7 yr _1 )(t Edd /10yr), and as a result, quiescent 
GCs, those with IMBHs now starved of fuel, should greatly outnumber active ones. 



Given a globular cluster space density of n gc ~ 4 Mpc 3 (Brodie & Strader 2006), 



we expect the density rate of flares with L = L Edd to be at most ~ 4000yr _1 Gpc~ 3 , if 
10 3 — 10 4 Mq black holes were prevalent in globular clusters. Therefore we would not yet 
expect to have detected such a flare. Observational constraints on the presence of IMBH 
would not then be stringent until a sufficiently large sample of galaxies should be sampled in 
order to reveal some of their GC members in a flaring state. Such objects should be searched 
for out to large distances. However, assuming L ~ 10 41 erg s -1 and T ~ T Edd , we find that 



EXIST (Grindlay 2004) would only be able to detect such a flaring event out to a distance 



of about 1 Mpc. 



4.2. ULX activity 

A further question is how fast the luminosity fades after the flare. This is important 
because we want to know whether the IMBH has faded below ULX levels before the next 
stellar disruption occurs. The answer to this question depends on how long it takes the last 
residues of the star to be ingested. One expects the infall rate M to decline as t" 5 / 3 for late 
times. Some material may, however, be stored for a longer time in an accretion disk. This is 
mainly because the specific angular momentum for Keplerian orbits grows oc r 1 / 2 , so angular 
momentum transport via disk viscosity requires that 10% of the debris goes out to 10 3 Rt 
and 1% to 10 5 Rt before being swallowed. 



The evolution of the bound debris onto a black hole for t > t a has been studied by 
Cannizzo et al. (1990) using a time dependent a-disk model. Cannizzo and collaborators 
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concluded that the black hole's luminosity fades more slowly: L oc i~ L2 and depends only 
weakly on the disk viscosity. The IMBH luminosity could then remain as high as Lx > 
10 39 erg s _1 for several hundreds of years after disruption. In a given object, this ULX 
activity would have a duty cycle of order 10~ 4 unless the rate of disruptions is much larger 
than 10~ 7 or the viscosity were so low that the bulk of the mass could be stored for many 
thousands of years at R ~ (where the dynamical timescale is only a few hours). 



4.3. Relevance to GCs and IMBH Growth 

The mass fraction that is ejected rather than ingested, although less spectacular than 
the accretion-powered flares could nonetheless have an important influence on the energy 
balance within a GC - similar than a supernova exploding in the same volume. When the 
star is disrupted in a single flyby, about half the debris is ejected in a fast moving spray of 
gas; the kinetic energy output being ~ 10 50 M^ g 3 erg. The ejecta would be braked and their 
kinetic energy thermalized, as they ran into the diffuse gas (probably originating from stellar 
mass loss). If this material is effectively shock heated, it could contribute to the X-ray and 
radio luminosity of the GC long after the star has been disrupted. Stellar disruption may 
have other implications - for instance, to the dynamics of line emitting regions within GCs 



(Chomiuk et al. 2008). 



The most tightly bound debris, on the other hand, would traverse an elliptical orbit 
before returning to R^. There would therefore be no conspicuous flare until, as discussed 
above, the bound debris fall back into the hole. Electron scattering opacity almost certainly 
dominates the radiative transfer and the photons will be trapped out to a radius R T = 
MK T /4:7rc. Energy will be dissipated at a supercritical rate as the material swirls closer to the 
hole, some gas may then be ejected in a radiation-driven wind. In principle it is possible for 
only a fraction e- l {RjR T ) ~ QMAe^Rj R & )- l (MjM & ) l / z Ml^ to be actually swallowed, 
all the remainder being ejected. If a fraction / of the mass from disrupted stars is accreted 
onto the black hole, the average rate of tidal disruptions required to form an IMBH out of 
a ~ 5OM progenitor would need to be > 10 _6 Mh i 3(//0.1)~ 1 (rGc/10 10 yr)~ 1 yr~ 1 , where tqc 
is the globular cluster age. In such cases, the density would, however, be high enough that 
runaway merging of high-mass main-sequence stars could lead directly to the formation of 



an IMBH (Baumgardt et al. 2006 Portegies Zwart et al. 2004). 
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Fig. 1. — A 1 M & solar-type star (modeled with more than 10 6 SPH particles) approaching 
a 10 3 M Q black hole on a parabolic orbit with pericenter distance -R m i n = Rt/3 is distorted, 
spun up during infall and then tidally disrupted. The panels show density cuts (in cgs units) 
through the orbital (xy-) plane before and after passage through pericenter. 
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specific energy [e] 



Fig. 2. — Differential mass distributions in specific energy for the 1 M stellar debris. At 
t =47.97 h, the amount of material unbound is ~ 47% of the initial mass of the star. 
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Fig. 3. — Density cuts (in cgs units) in the orbital plane of the bound debris 22.9 hours after 
disruption. The most tightly bound debris would transverse an elliptical orbit with major 
axis ~ 10 4 i? g before returning to R ~ R p , where radial focusing of orbits acts as an effective 
nozzle. These orbits are focused back into the original orbital plane at pericenter. This 
causes the formation of a pancake shock that weakly redistributes the orbital parameters 
and damps out some of the vertical motion. 
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Fig. 4. — The rates at which the stellar debris returns to the vicinity of the black hole for 
different types of stars: A [1,1,0.6] M & solar-type star modeled with T = [5/3, 1.4, 1.4] shown 
as a solid [blue, red, black] line. The infalling matter at a rate that drops off roughly as t~ 5 / 3 
(gray line). Inset: the internal density structure of the various stars. 



